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Abstract 

Effect of phase fluctuations on superconducting states with anisotropic order 
parameters is studied in a BCS like lattice model of cuprate superconductors. The 
degradation of the mean field transition temperature due to phase fluctuations is 
estimated within a Kosterlitz-Thouless scenario. Values of the interaction param- 
eters for optimal doping, corresponding to a stable superconducting state of S xy 
symmetry, which fit the nodal structure of the superconducting order parameter 
in the Bi2212 compound, are obtained. The angular position of the node is found 
to be insensitive to the dopant concentration. 
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The symmetry of the order parameter (OP) in high temperature superconductors 
continues to be a matter of much controversy. A series of experiments [1-6] have indicated 
that the OP in high temperature superconductors is highly anisotropic J7|. However, 
these experiments have conflicting features [8,9] which have not yet made it possible 
to determine the symmetry of the OP in an unambiguous way. Recent angle resolved 
photoemission spectroscopy (ARPES) experiments 0, |6| have been able to shed new 
light on the symmetry of the OP, owing to their high resolution (~10 meV) compared 
to the size of the superconducting gaps in the cuprates (~50 meV). For high quality 
Bi2212 single-crystal samples, the gap on the Fermi surface (FS) vanishes at two points 
symmetrically displaced from the 45° direction in the Brillouin zone by about 10 degrees 
U at the optimum dopant concentration (5 = 0.17) which has the highest Tc. Such 
a nodal structure is inconsistent with the d x 2_ y 2 symmetry, which has only one node 
per quadrant along the 45° direction. Also, a small admixture of an isotropic S-wave 
component to the d x 2_ y 2 OP does not increase the number of nodes but merely displaces 
the node away from the 45° direction. 



Norman and coworkers have, in a recent set of papers |1(J , made a strong case for an 
OP whose symmetry is predominantly S xy . They considered a phenomenological model 
of electrons on a square lattice with an on-site repulsion (Vo) and attractive nearest 
neighbour (Vi) as well as next nearest neighbour (V2) interactions. The Hamiltonian of 
the system is given by H = H + Hi where 



and 



1.(7.(7 
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The band dispersion was chosen to fit the normal state quasiparticle dispersion ob- 
tained from the ARPES data on Bi2212 single crystals [|K|]. A six parameter tight 
binding Hamiltonian with [t ,...,t 5 ) = [0.131,-0.149,0.041,-0.013,-0.014,0.013] (in eV) 
was found to reproduce the observed quasiparticle dispersion. Here t is the orbital 
energy, t\ nearest neighbour (nn), t 2 next nearest neighbour (nnn) etc. hopping matrix 
elements. 

A mean field analysis of the spin-singlet pairing instabilities of the above model 
pTOfl showed that the strongest instabilities correspond to OP's with S xy and d x 2_ y 2 
symmetries because they best exploit the large single particle density of states (DOS) 
just below the FS; their relative stability is determined by the ratio of V\ and V%. A 
consistent fit to the ARPES data could be obtained in a range of values of V\ and V2 
while Vo was constrained to be small. 

However, the above calculations overlook the strong fluctuations that exist in the 
cuprates owing to their quasi two dimensional nature. In this letter, we study the phase 
diagram of the Hamiltonian in Eq.(l) with the incorporation of the effects of fluctuations 
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in the phase of the OP. We calculate the superfluid phase stiffness or helicity modulus 
(p s ) below the mean field transition temperature (T C MF ), and use it in conjunction with 
the Kosterlitz-Thouless (KT) relation p s (T^ T ) = \k B T^ T to determine the reduced 
transition temperature (T^ T ). The interplanar Josephson coupling, that exists in the 
cuprates, does not lead to a strict KT like scenario as is being assumed here (except in 
thin films). However, the magnitude of T^ T is determined by p s and can therefore be 
expected to yield a good estimate of T c degradation due to phase fluctuations. With this 
improved bound on T c , we determine the relative stability of superconducting OP's with 
Bi (d x 2_ y 2), B 2 (d xy ) and Ai (a 3 dimensional representation comprising S, S* and S xy ) 
symmetries. We also study the nodal structure of the gap function on the FS. 

We find that the region of stability of the Ai (S xy ) OP with respect to d x 2_ y 2 and 
d xy , is reduced compared to the mean field result. An OP which is predominantly S xy 
with very small isotropic S-wave (S) and extended S-wave (S*) components correctly 
reproduces the gap structure predicted from ARPES measurements. The best fit, cor- 
responding to a T* T ~ 100K and a stable Ai phase with a nodal structure consistent 
with ARPES measurements, occurs for V < 15meV, V\ in the range 30 — A5meV and 
Vi around 70 — 75meV. For these parameters we find the variation of the transition 
temperature (T^ T ) to be qualitatively correct while the angular positions of the nodes 
on the FS are nearly independent of 8. 

Standard mean field factorization of the Hamiltonian (Eq.(l)) gives the gap equation 
Ar = - V V(k - k') tanh ( (2) 

where the quasiparticle energy is E% = yj (eg — e^) 2 + | Ag| 2 and Ag is the BCS gap 
function. The pairing interaction, V(q) — Vq + 4Vi(cosq x + cosq y ) + 8V2 cos q x cos q y , is 
separable and can be written as V(k — k') = X) i=0 ViViWViik') where r) (k) = 1, rji(k) = 
|(cos k x + cos k y ), r)2(k) = |(cos k x — cos k y ), i] 3 (k) = cos k x cos k y , rj^k) = sin k x sin k y , 
corresponding to [S, S*, d x 2_ y 2, S xy , d xy ] symmetries, and (V0...V4) = (Vo, 8V1, 8V1, 8V2, 
8V2). We have ignored terms corresponding to triplet pairing. It is easy to see that 

— # 

upon writing the gap functions as Ag = J2i Vi(k)Ai, the linearized gap equation (whose 
non-vanishing solution first appears at T C MF ) factorizes into two independent equations 
for A2 and A4 (Bi and B2 representations) and three coupled linear equations involving 
A , Ai, A 3 (Ai representation). Since A 3 is the largest component within the Ai repre- 
sentation, the OP symmetry is predominantly S xy . We ignore the possibility of further 
transitions at lower temperatures. 

We now proceed to the calculation of the p s and T^ T . We first solve the gap equations 
below T C MF , choosing the gap function Ag = J2ien \Vi(k) ■ As an illustration, we plot 
the gap functions A , Ai, A 3 (Ai) and A 2 , A 4 (Bi and B 2 respectively) in the inset 
of Fig.l, as a function of temperature for optimal parameters (see below). To calculate 
p s , we introduce a transverse vector potential and choose the gauge A y = 0. Then, the 
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hopping matrix elements in Hq (Eq.(la)) acquire an additional phase factor according 
to the Peierls substitution t i3 - — > exp[^ f2* A ■ dl). 



The electron current operator j x (Ri) can be found (to linear order in A x ) by differ- 
entiating Hq with respect to A x (Ri) and is of the form 

dH 



j x (Ri) = -c- 



f x ara (Ri 



+ Jx 



dia {Ri 



(3) 



dA x (Ri 

The first term on the RHS of Eq.(3) is of zeroth order in A x and is the electron velocity 
operator. It represents the paramagnetic response of the quasi-particles. The second 
term is of linear order in A x and its physical origin is related to the Meissner screening 
of the condensate. The average diamagnetic current is obtained to be 



■dia 
Jx 



(q) 



._^__LW c t c , 

k.a 



AM 



(4) 



where the angular brackets correspond to an average in the mean field superconducting 
state. From Eq.(4) we see that the diamagnetic contribution to the phase stiffness (pf a ) 

is proportional to j^Z^^cl c^Jj^r, the mean electronic kinetic energy along the 
x-direction JLT|. Unlike models in which the electron dispersion is parabolic, pf a for a 
lattice superconductor is not proportional to the average electronic density. Thus the 



use 



12[| of the continuum expression pf a oc — J2% a 



9e k\ gKf£-!£) 



(where / is the Fermi 



function) is incorrect and could leave out important physics. 

From linear response theory, the long wavelength paramagnetic current is found to 



be 



f x ara (q) 



lim lim K xx (q,uj) 



A x {q) 



(5a) 



where K xx (q, to) — — % J dt6(t)e lU)t ([j x ara (q,t), j x ara (—Q,0)])- We work in a transverse 
gauge and vertex corrections required to get a gauge invariant current have been ignored. 
The correlation function in Eq.(5a) is easily evaluated and found to be 

2 



lim lim K xx 



(2» = ^E(^ 



de % Ydf{E- k 



dk x 



dE t 



(56) 



From Eqs (4) and (5) we obtain the superfluid phase stiffness 

2 



Ps{T) 



1 

2iV 



E 




dm) 



+ 



2 dk 2 x 



1 



E, 



tanh 




(6) 



In Fig.l we show p s as a function of temperature for various OP symmetries. The 
results indicate that the temperature dependence of p s does not show any major qual- 
itative change for the various types of OP's considered. Thus the fluctuation corrected 
phase boundaries are not expected to be very different from the mean field case. 
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The KT transition temperature (T c ) is found by solving 

Ps(T? T ) = ~T C KT (7) 

7T 

We use the expression for p s (Eq.(6) in Eq.(7)), assuming mean field values of the OP, 
and obtain T^ T . This procedure was employed earlier by Denteneer et. al. Jl3 in a 
study of the 2-d attractive Hubbard model and gives a qualitatively correct picture even 
in the limit of strong coupling where T^ T is inversely proportional to the coupling. Since 
we neglect the renormalization of the OP (and hence p s ), due to the presence of vortex 
like fluctuations, the discontinous drop in p s , characteristic of KT transitions, is not 
obtainable here. Hence our results give an upper bound on the true KT transition tem- 
perature. Another possible difficulty is that the expression in Eq.(7) has been derived 
for isotropic S-wave superconductors. However, even in the anisotropic case, the asymp- 
totic form of the vortex-vortex interaction (at high vortex density) is logarithmic |14| 
and so corrections to Eq.(7) are likely to be small. A comparison of T C XT thus obtained 
for different OP's enable us to determine the phase diagram of the model Hamiltonian. 

We now discuss our results. For the sake of clarity we present the phase diagrams 
in the Vq — V 2 and V\ — V 2 planes only. In Fig. 2 we have drawn the phase diagram 
in the V$ — V 2 plane (we set V\ = 0) for several values of 5. For comparison, we also 
include the mean field phase boundaries [fLOfl . The dopant concentration (5) indicates 
the deviation from half-filling. It is clear from Fig. 2 that fluctuation effects are negligible 
at weak coupling (|T^| < 50meV). To obtain a stable S xy (Ai) solution at low filling, 
V is required to be rather small (V < 0.4 eV at 5 = 0.17), even smaller than the value 
indicated by the mean field analysis and larger values of | V2I/V0 are required to stabilize 
the S xy state. At high values of 5, the Ai phase is stable for much larger values of V (a 
few eV) provided [V^l < 90meV. This is understandable as the onsite repulsion, which 
disfavours the Ai state, has the largest effect close to half filling and becomes relatively 
ineffective at low filling. 

In Fig. 3 we plot the phase boundary separating the d x 2_ y 2 and Ai superconducting 
states in the V\ — V 2 plane (setting V = 0). Fluctuation effects are once again seen to 
be negligible at weak coupling but lead to a larger critical value of V 2 /Vi (to stabilize 
the S xy state) at intermediate and strong couplings. Both the mean field and KT 
phase boundaries are insensitive to 5. This is not surprising as both the d x 2_ y 2 and S xy 
solutions have a very similar pairing density of states JjlJ and thus the phase boundary 



is determined almost entirely by the interaction strength. 

We have found the values of interaction parameters at 5 = 0.17 (optimal doping) 
that produce a stable superconducting state with the Ai symmetry and a transition 
temperature T^ T ~ 100 K whose OP has the observed nodal structure. We find that 
the on-site repulsion Vq is constrained to be small (< l5meV), V 2 lies in the range 
70 < —V 2 < 75meV whereas V\ can be chosen in a relatively broader range 30 < 
—V\ < 45meV. In Fig. 4 we have plotted the angular dependence of |At(T = 0)| on 
the FS (0 measures the angular deviation from the line passing through the (0, 7r) and 
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(n, 7r) points), for typical values of Vo, V\ and V%, lying within this range. A puzzling 
feature of these results, which does not have any apparent microscopic justification, is 
the extremely small value of Vo, whereas the actual bare Cu onsite repulsion is known 
to be about a few eV. 

We have studied the variation of T^ T with 5 for the optimal parameters determined 
above. We find that T^ T has the correct qualitative behaviour (see inset of Fig. 5) with 
the highest Tjf T being achieved for 5 ~ 0.25. This is not surprising as the pairing 
density of states (for a superconducting instability in the Ai state) at the Fermi energy 
is peaked close to this value of 5 flO|| . We have also found the nodal structure of the 



OP as a function of 5. The magnitude of the gap function on the FS away from the 
nodes is strongly dependent on 5 (reflecting the variation in Tjf T ). However the angular 
positions of the nodes on the FS are insensitive to 5 (see Fig. 5 and its inset) and stay 
close to their positions for 5 = 0.17. 

In conclusion, we have studied the competition between superconducting OP's with 
different symmetries within a phenomenological BCS model and obtained a phase dia- 
gram which incorporates the effects of fluctuations. Our phase diagram, together with 
the magnitude of the transition temperature and the nodal structure of the OP, enable 
us to identify the best choice of parameters appropriate to the optimally doped Bi2212 
compound. We find that the transition temperature has the correct doping dependence 
and the angular positions of the nodes of the gap function on the FS are nearly inde- 
pendent of S for the best parameters. 
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FIGURE CAPTIONS 



Fig.l. p s as a function of temperature at 8 = 0.17, Vq = 15meV, V\ = —AOmeV and 
V2 — —75meV for the A 1; B x and B 2 symmetries. The intersection of the straight line 
with p s gives T C XT . Inset shows the variation of the mean field OP's with temperature. 

Fig. 2. The mean field (dashed lines) and KT (solid lines) phase boundaries indicating 
the regions of stability of the Ai and d xy phases are shown in the (Vo, — V2) plane for 
various 8. The inset shows the KT phase diagram for 8 — 0.17 on an expanded scale. 

Fig. 3. The mean field and KT phase boundaries indicating the regions of stability of the 
Ai and d x 2_ y 2 phases are shown in the (— V2, — VI) plane for various 8. The inset shows 
the KT phase diagram for 8 = 0.17 on an expanded scale. 

Fig. 4. The angular dependence of the gap function |Ak(T = 0)| on the Fermi surface is 
shown for 8 = 0.17 and various values of interaction parameters in the optimal range. 
The angles are measured with respect to the (n, ir) point as in ref. [10]. 

Fig. 5. |A k (T = 0)| on the Fermi surface with V = 15meV , V\ = —40meV and 
V 2 = —74meV (optimal parameters) for various 8. Curves 1-5 correspond to 8= 
0.10,0.13,0.17,0.21 and 0.26 respectively. Inset shows the angular position (a) of the 
node on the Fermi surface (measured from the 45° direction) and T^ T (in °K) as a 
function of 8 for optimal parameters. 
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